Centrality and pseudorapidity dependence of charged hadron production at 
intermediate px in Au + An collisions at ^/s^ = 130 GeV 
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We present STAR measurements of charged hadron production as a function of centrality in Au + 
Au collisions at ^s„„ = 130 GeV. The measurements cover a phase space region of 0.2 < pr < 6.0 
GeV/c in transverse momentum and — 1 < ?7 < 1 in pseudorapidity. Inclusive transverse momentum 
distributions of charged hadrons in the pseudorapidity region 0.5 < < 1 are reported and 
compared to our previously published results for \ri\ < 0.5. No significant difference is seen for 
inclusive pr distributions of charged hadrons in these two pseudorapidity bins. We measured dN/dr} 
distributions and truncated mean pr in a region of pr > Pt'*, and studied the results in the 
framework of participant and binary scaling. No clear evidence is observed for participant scaling of 
charged hadron yield in the measured pr region. The relative importance of hard scattering process 
is investigated through binary scaling fraction of particle production. 

PACS numbers: 25.75.Dw, 25.75.-q 
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I. INTRODUCTION 



pseudorapidity due to change of momentum is 



Quantum Chromodynamics (QCD) is considered to be 
the underlying theory of the strong interaction which 
governs hadron production in nuclear collisions. The 
strong interaction is usually divided into soft processes, 
which involve small momentum transfer, and hard pro- 
cesses, which can be calculated using perturbative QCD. 
The Relativistic Heavy Ion Collider (RHIC) experiments 
at the Brookhaven National Laboratory investigate prop- 
erties and evolution of matter at high temperature and 
energy density. At RHIC energies, the hard processes be- 
come more evident in comparison to previous heavy ion 
experiments and can be used to probe the early state of 
the collision system. A high energy parton produced via 
hard scattering may lose energy in the hot /dense medium 
through gluon bremsstrahlung and multiple scatterings 
before hadronization P, Q , leading to a suppression of 
high Pt hadron production. The magnitude of the en- 
ergy loss provides an indirect signature of QGP forma- 
tion. Since parton energy loss is directly proportional to 
gluon density, the energy loss would be much larger in a 
partonic medium than in hadronic matter Q. 

Partonic energy loss can be investigated through com- 
parison of hadron yield as a function of pr between 
nucleus- nucleus collisions and p -\- p or p + p collisions. 
In order to do so, scaling factors which account for the 
nuclear geometry, the number of participant nucleons, 
Npart, and the number of binary nucleon-nucleon colli- 
sions, Njjim are calculated from theoretical models. Ex- 
perimental results from RHIC, including our earlier anal- 
yses in the pseudorapidity region I77I < 0.5, have indi- 
cated a suppression of hadron production for pT > 2 
GeV/c in central Au -I- Au collisions relative io p+p and 
p+p collisions j^lHI^' This is in contrast to the SPS re- 
sult from central Pb -I- Pb collisions at ^s„„ — 17 GeV, 
which shows an excess of i:^ production for 2 < < 4 
GeV/c 0,11. The RHIC measurements are striking con- 
sidering that known nuclear effects, like the Cronin ef- 
fect U and radial flow ts'^d to enhance hadron 
yields at high px- The RHIC results for high pT hadron 
suppression agree qualitatively with calculations based 
on fragmentation models, which attribute the high pT 
hadron suppression to medium induced parton energy 
loss 11]. 

Another known nuclear effect, nuclear shadowing, 
also modifies particle production at high px- Calcula- 
tions of this effe ct [ia | based on the EKS98 shadowing 
parametrization [l3j predicted it to be small in the pt 
and pseudorapidity region covered in this measurement. 
However, another study 0| found a much larger shad- 
owing effect for heavy nuclei at RHIC. Therefore, a mea- 
surement of particle production as a function of pT and 
pseudorapidity may provide a constraint on the shadow- 
ing effect. 

Partonic energy loss may also be studied by the pseu- 
dorapidity dependence of hadron production. Change of 
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The pseudorapidity distributions would be modified as a 
result of the parton energy loss if the momentum change 
rate {5p/p) due to the energy loss is different along the 
transverse and longitudinal direction. In addition, Pol- 
leri and Yuan (l^] pointed out that the degree of the en- 
ergy loss may also depend on the pseudorapidity region 
in which a jet is produced because the energy loss is pro- 
portional to the particle density in pseudorapidity. The 
pseudorapidity dependence of high pT hadron production 
provides a means to probe the initial density of matter 
along both the transverse and longitudinal directions. 

In this article, we present measurements of hadron pro- 
duction in Au -I- Au collisions at yfs^ = 130 GeV as a 
function of centrality, p^ and ry. In Sec. II we will briefly 
describe the STAR experimental setup and then give a 
description of data analysis techniques that were used to 
obtain the inclusive transverse momentum distributions 
for charged hadrons. We will also discuss the parameter- 
ization of inclusive transverse momentum distributions 
in p + p collisions at -/s = 130 GeV and the calcula- 
tions of Npart and Nbin- In Sec. HI results from the data 
analysis will be reported and compared with model cal- 
culations. Physics implications of our measurements are 
discussed in Sec. IV. And we will then summarize our 
measurements in Sec. V. 



II. ANALYSIS 
A. Experimental Setup and Data 

Measurements presented in this article are based on 
two data sets of Au + Au colhsions at y/s^ = 130 GeV, 
which were recorded by the STAR detector at RHIC. A 
detailed description of the STAR detector can be found 
elsewhere '16^ . The two data sets comprise minimum bias 
and central collision triggered events which correspond 
to approximately the most central 10% of the Au + Au 
geometric cross section. Charged particle tracks of an 
event were detected in the Time Projection Chamber 
(TPC) with a pseudorapidity coverage |?7| < 1.8 and com- 
plete azimuthal symmetry. The transverse momentum 
of a track is determined by fitting a circle through the 
transverse coordinates of the primary event vertex and 
the space points along the track in the TPC. The total 
momentum can be calculated using this radius of curva- 
ture in a 0.25 T magnetic field and the polar angle of the 
track. The procedure involves a three dimensional fit us- 
ing three coordinates of the primary vertex determined 
from all of the tracks reconstructed in the TPC. The 
primary vertex position along the beam direction, Zytx, 
has a wide spread with one standard deviation about 100 
cm. To increase detection efficiency of the tracks within 
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\r]\ < 1, we required the events to have a primary ver- 
tex \zytx\ < 75 cm. After the event selection cuts, the 
minimum bias data set contained ~ 181k events and the 
central data set contained ~ 365k events. 

Centrality selection is based on the uncorrected pri- 
mary charged particle multiplicity Net within |?7| < 0.75 
and pt < 1.5 GeV/c. The requirement on 77 range max- 
imizes the number of tracks used to define centrality in 
an event while keeping the tracking acceptance approxi- 
mately constant. The percentage of the geometric cross 
section is determined in the same way as that published 
by STAR previously 0, where the negatively charged 
hadron multiplicity Nf^- distribution within \r]\ < 0.5 
was used. The data set is divided into seven centrality 
bins, and the most central bin is 0-5% (the top 5% of the 
multiplicity distribution) while the most peripheral bin 
is 60-80%. 

The analysis in this article covers a transverse momen- 
tum region of 0.2 < px < 6.0 GeV/c. Accepted primary 
tracks have \ri\ < 1, at least 25 space points in the TPC 
used in the track fit out of 45 pad rows, a fit probability 
of being a primary track greater than 0.05, and a distance 
of closest approach to the primary vertex less than 1 cm. 
These track quality cuts were varied to estimate the sys- 
tematic uncertainty. Acceptance and efficiency were de- 
termined by embedding simulated tracks into actual Au 
+ Au collision events. 

The measured high px hadron yield is sensitive to small 
spatial distortions of the TPC alignments in both az- 
imuthal and longitudinal directions. A measurement of 
the summed hadron yield, {h~^ + h~)/2, is less sensitive to 
such distortions than the yield of one charge sign alone. 
Using 12 sectors from each of the TPC ends as indepen- 
dent detectors for high px hadrons, we estimated the sec- 
torwise (azimuthal direction) variations of the yields to 
be less than 5%. The variation of the yield between the 
hadrons crossing and not crossing the central membrane 
of the TPC was found to be approximately proportional 
to Pt with a value of 11% at px = 5.5 GeV/c. The typi- 
cal correction factors for the acceptance and efficiency are 
given in Table^as "Tracking" . The systematic uncertain- 
ties incorporate acceptance, efficiency, track quality cuts, 
and the effects of the spatial nonuniformity. The tracking 
and other correction factors and their systematic uncer- 
tainties given in Table Q for I77I < 0.5 differ from those 
given in our previous paper 4] because different track 
quality cuts and other correction procedures were used. 

Finite momentum resolution tends to spread particles 
to neighboring bins in a momentum histogram, especially 
for an exponentially falling spectrum. This smearing ef- 
fect cannot be neglected at higher px where the momen- 
tum resolution is limited by the strength of the magnetic 
field and the TPC spatial resolution. We used the em- 
bedding technique to determine the pr resolution. For 
Pt > 0.5 GeV/c within \r]\ < 0.5 the Gaussian distribu- 
tion of track curvature fc cx 1/pt has a relative width of 
Sk/k = 0.013 -I- 0.015pT/(GeV/c) for central events and 
5k/k = 0.012 + 0.012pT/(GeV/c) for peripheral events. 



Within 0.5 < |?7| < 1, 5k/k = 0.014 0.010pT/(GeV/c) 
for central events and 5k/k = 0.014 0.0072]5T/(GeV/c) 
for peripheral events. 

The fact that the pT resolution for 0.5 < It^I < 1 is 
better than that for I77I < 0.5 is due to the competition 
between two opposing effects. For a given pT track in 
the TPC, the hadron with higher 77 tends to have fewer 
space points, hence poorer resolution, but shorter drift 
distance, hence better resolution. 

The magnitude of the pt resolution determined from 
the embedding technique did not include the effect of the 
primary vertex resolution. The effects of the pt smearing 
due to the primary vertex resolution, to the charge-sign- 
dependent distortion, and to the weak decay background 
tracks, have been empirically derived from the compar- 
ison between real and embedded tracks. The combined 
effect within 0.5 < < 1 was found to be larger than 
that within \ri\ < 0.5. This is partially due to the fact 
that the magnitude of the charge-sign-dependent distor- 
tion in the higher i] region is larger. 

The two contributions have been convoluted into a 
power law function to fit the data, and then the ratio 
of the fitted function to its convoluted one gives the pt 
smearing correction factor j^^. Because the two con- 
tributions have opposite \r]\ dependence, the overall pt 
smearing correction factors for the two 77 regions happen 
to be comparable. The typical pt smearing correction 
factors and their systematic uncertainties are also given 
in Tabled 



B. Background 

The most significant backgrounds for the high pT 
charged hadron yield as seen in Tabled come from parti- 
cle weak decays and antinucleon annihilation in detector 
material. The contamination rate for each background 
source was estimated using detector response simulations 
with events generated by the HIJING model How- 
ever, the Pt dependence of production of weakly decay- 
ing particles, primarily Kg, A, A, and of antinucleons, p, 
n, in HIJING is not consistent with experimental mea- 
surements. We corrected those predicted yields usin g th e 
measured spectra of p 0, 113, A and A [23 , and JtT" |2^ , 
together with those of h~ |l8ll2ll|. for pT < 2.4 GeV/c 
in the mid-rapidity region in the most central bin. The 
corrections used in calculating the background fractions 
are shown in the upper panel of Fig. The curves are 
polynomial fits to the data points and are used in the 
interpolation due to different px binning. For pT > 2.4 
GeV/c we simply assumed the yield ratios to be constant. 
Systematic uncertainties of 50% and 100% of the over- 
all background fraction are assigned for the regions of 
Pt < 2.4 GeV/c and px > 2A GeV/c, respectively p^ . 

The contamination rate for all background sources 
shows almost no centrality dependence from the Monte 
Carlo HIJING events. Therefore, centrality dependence 
of the background fraction is mainly determined by the 
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FIG. 2: 1] acceptance correction function from model calcu- 
lations: ratios of pr spectra in two different rj regions to that 
within \rj\ < 2.5, in which the UAl Collaboration published 
its inclusive charged particle pr spectra. 
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FIG. 1: Measurements used in background studies. Up- 
per panel: ratios of the measured pT yield ratios to those 
of HIJING in the most central bin. Lower panel: measured 
rriT inverse slope parameters as functions of centrality rep- 
resented by measured negatively charged hadron multiplicity 
Nf^- within \ri\ < 0.5. Curves are polynomial fits to data 
points. 



measured spectra in various centrality bins. In the lower 
panel of Fig. ^ we show the measured transverse mass 
{rriT = \Jv\ + "^o) inverse slope parameters of exponen- 



tial fits to p [23, A and A |23j, and K'^g |2J| spectra in 
the mid-rapidity region as functions of the measured neg- 
atively charged hadron multiplicity T^h- within \'r]\ < 0.5. 
We use these to correct for different centrality binning in 
our analysis. The polynomial fits are used to interpolate 
the inverse slope parameters in the centrality bins used 
in this analysis. 

Pseudorapidity dependence of the background fraction 
is studied using the Monte Carlo HIJING events. For 
Pt < 2 GeV/c the ?7-dependence of backgrounds is neg- 
ligible within —1 < 77 < 1 while for pr > 2 GeV/c the 
background fraction decreases with pT and \r]\. For exam- 
ple, at Pt = 5.5 GeV/c the background fraction predicted 
within 0.5 < |ry| < 1 is only 40% of that within |?7| < 0.5. 
The typical background correction factors and their sys- 
tematic uncertainties are given in Table |l| The total 
systematic uncertainties of the measured spectra within 
I77I < 0.5 (0.5 < \r]\ < 1) at the highest bin pT = 5.5 
GeV/c are « 24% (« 18%) for central events and « 17% 
15%) for peripheral events. 



In the absence of any NN collision data at ^/s — 
130 GeV, a NN reference spectrum is obtained by ex- 
trapolation of the UAl p + p data for ^/s = 200 — 
900 GeV m. The UAl inclusive charged particle pt 
spectra within \r]\ < 2.5 were fitted by the pQCD in- 
spired power law function 
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The fit parameters were used to extrapolate to our 
energy, giving Co-j„ = 267tg mb/(GeV/c)2 {ain de- 
notes the inelastic cross section of NN collisions), 
po = l-90^o;o9 GeV/c, and n = l2Mtof^ at = 
130 GeV 3. The superscripts and subscripts are curves 
that bound the systematic uncertainty. 

However, the UAl acceptance is different from 
star's. Corrections were made to the UAl reference 
for our 77 acceptance based on two independent pQCD 
calculations: PYTHIA 26] and Vitev 123 • When the 
K factor in PYTHIA is set to 1.5, PYTHIA calculations 
for 200 GeV p + p collisions are in reasonable agreement 
with the UAl measurement of the inclusive charged par- 
ticle Pt spectrum and with the UA5 measurement 
of the pseudorapidity density distribution Similar 
PYTHIA calculations are in reasonable agreement with 
the STAR measurement of the inclusive charged hadron 
Pt spectrum within \rj\ < 0.5 for p + p collisions at 
y/s = 200 GeV . Fig. El shows the pT-dependent cor- 
rection functions for two 77 regions at ^/s = 130 GeV, 
obtained by averaging over the two pQCD calculations. 
The solid curve is the ratio of d'^N / dpTdrj within jr/j < 0.5 
to that within jryj < 2.5, and the shaded area shows its 
systematic uncertainty. The dot-dashed curve shows the 
same ratio for 0.5 < jryj < 1, and the similar magnitude of 
the uncertainty on the ratio of 0.5 < jr;! < 1 to \ri\ < 2.5 
is not shown. Multiplicative corrections of 1.35 ± 0.09 
and 1.33 ± 0.09 at pt — 5.5 GeV/c have been obtained 
for I77I < 0.5 and for 0.5 < \ri\ < 1, respectively. The 
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TABLE I: Typical multiplicative correction factors and systematic uncertainties, applied to the yields for peripheral and central 
collisions within I?;] < 0.5 and within 0.5 < < 1. 



Pseudorapidity 




PT = 2 GeV/c 


PT = 5.5 GeV/c 


Centrality 


60-80% 


0-5% 


60-80% 


0-5% 


|r/| < 0.5 


Tracking 
Pt Smearing 
Background 


1.16 ±0.10 
1.01 ±0.01 
0.92 ± 0.04 


1.71 ±0.15 
1.00 ±0.01 
0.88 ±0.06 


1.22 ±0.16 
0.89 ± 0.02 
0.90 ±0.10 


1.65 ±0.22 
0.70 ± 0.06 
0.85 ±0.15 


0.5 < < 1 


Tracking 
Pt Smearing 
Background 


1.29 ±0.11 
1.01 ±0.01 
0.92 ± 0.04 


1.78 ±0.15 
1.01 ±0.01 
0.88 ± 0.06 


1.31 ±0.18 
0.89 ± 0.02 
0.96 ± 0.04 


1.71 ±0.23 
0.72 ±0.07 
0.94 ±0.06 



difference between jryj < 0.5 and 0.5 < jTyj < 1 is quite 
small, indicating a relatively fiat rj distribution within 
— 1 < ?7 < 1 for a broad px region. The STAR measure- 
ment [3 is consistent with the UAl p + p data for 200 
GeV after applying a similar 77 acceptance correction. 

We derived ain in the NN reference at y/s = 130 GeV 
of 40 ± 3 mb by requiring dN/dridrjl < 0.5), which was 
obtained by integrating the extrapolated spectrum af- 
ter applying the rj acceptance correction, to be 2.25, 
which was determined from the energy dependence of 
dN/dr]{Tj = 0) J9]. 



D. Participant and Binary Collision Determination 

The number of participant nucleons, Npart, and the 
number of binary nucleon-nucleon collisions, Nun, in a 
nucleus-nucleus collision are used to compare experimen- 
tal results with model predictions. Unfortunately, at 
RHIC Npart and N^in cannot be measured directly and 
have to be obtained in a model-dependent way. Consid- 
erable discrepancy exists among various model calcula- 
tions, especially for peripheral collisions ^30|. 

We first investigate Npart and obtained from a 

Monte Carlo (MC) Glauber model calculation j4,, 31] . In 
the Monte Carlo Glauber model, each of the nucleons in a 
nucleus A is randomly distributed using a Woods-Saxon 
nuclear density distribution 

= T, rP' — wra ' (2) 

1 ± exp[(r - ro)/D\ 

with normalization to J p{r)dr = A and parameters: nu- 
clear radius rg and surface diffuseness D. All nucleons in 
either nucleus for a nucleus-nucleus collision are required 
to be separated by a minimum distance. The calculated 
da / dNpart or da /dNun distribution was divided into bins 
corresponding to common fractions of the total geometric 
cross section to extract the average Npart or Nbin for each 
centrality bin. The systematic uncertainties on Npart and 
Nhin were estimated by varying the Woods-Saxon param- 
eters, by varying the ain value, and by including a 5% 
uncertainty in the determination of the total geometric 
cross section. 



We also investigate calculations of Npart and Nun us- 
ing two dynamic models, HIJING ^ and VENUS [13 . 
We compare these calculations with results from the 
Monte Carlo Glauber model calculation to shed light on 
the model dependent uncertainties of Npart and Nhin- 

The VENUS model is based on the Gribov-Regge the- 
ory and string fragmentation. The HIJING generator is 
an example of a two-component model: the momentum 
transfer of the soft process is treated phcnomenologically 
and the hard processes are calculated by pQCD. The ex- 
cited nucleons after collisions are stretched out as quark- 
diquark strings and fragment based on the Lund frag- 
mentation scheme '3^^. The parton energy loss in dense 
medium (quenching) and nuclear modification of parton 
structure functions (shadowing) are also modeled in HI- 
JING. 

Both dynamic models describe nuclear collision geom- 
etry using the Woods-Saxon nuclear density distribu- 
tion and the eikonal formalism to determine the prob- 
ability for each binary nucleon-nucleon collision, and to 
compute Npart and Ntin ■ Table |n] shows the compar- 
isons of the nuclear geometries implemented in HIJING, 
VENUS, and the Monte Carlo Glauber model for Au + 
Au collisions at y/s„„ = 130 GeV. The overlap func- 
tion, which defines the probability for a nucleon-nucleon 
collision at a given impact parameter 6, has the form of 
l-exp(-2r2(6)) in HIJING with n{b) defined in Table|TT| 
{po — 3.9 and 7769(5) = asoft{s) /2) while it is a step 
function, e{R - 6), in VENUS and MC Glauber. 

The correspondence between the centrality classes de- 
fined by measured charged particle multiplicity and those 
defined by modeled impact parameter was used to extract 
the average Npart and Nun from these dynamic models 
for a given centrality bin. Variations of average Npart and 
Nbin for different centrality selections were estimated us- 
ing the Monte Carlo events from the HIJING model. The 
event classes corresponding to the same fractional cross 
section were selected by cuts on b, Nch , Npart , and N^in ■ 
The average Npart and N^in by different cuts in HIJING 
are consistent within 2% for each centrality bin except 
the 60-80% most peripheral bin, where the discrepancy 
is at a level of 6%. 

The results of Npart and N^in from the models are 
shown in Fig. O as their ratios to Npart or Nbin from 
the Monte Carlo Glauber calculation. The participant 
scaling exponents a, which are defined in the expression 
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TABLE II: Comparisons of nuclear geometries implemented in various models for 130 GeV Au + An collisions. 



Model 


HIJING 1.35 


VENUS 4.12 


MC Glauber 


Woods-Saxon 
Parameters 


ro = 6.38 fm 
D = 0.535 fm 


ro — 6.64 fm 
D = 0.540 fm 


ro = 6.5 ±0.1 fm 
D = 0.535 ± 0.027 fm 


Minimum Distance 
of Two Nucleons 


0.4 fm 


0.8 fm 


0.4 fm 


Nucleon-Nucleon 
Overlap Function 


n{b) = (1 + (7jet/asoft)X0i0 


9{R-b) 


e{R - b) 


Maximum Impact 
Parameter 


25.6 fm 


24.1 fm 


no restriction 


Nucleon-Nucleon 
Cross Section ain 


38.7 mb 


37.4 mb 


41 ± 1 mb 


Total Geometric 
Cross Section 


7.27 b 


7.34 b 


6.9 ±0.4 b 



• HIJING /Glauber 
A VENUS /Glauber 
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FIG. 3: Ratio of the number of participants Npart (upper 
panel) or the number of binary collisions Ntin (lower panel) 
determined from different models to that from a Monte Carlo 
Glauber calculation. Shaded areas show the uncertainties 
of Npart or Nbin from the Monte Carlo Glauber calculation. 
Curves are to guide the eye. 



FIG. 4: Inclusive pr distributions of (ft^ )/2 within 0.5 < 
\ri\ < 1. Noncentral bins are scaled down by the indicated 
factors. The combined statistical and systematic errors are 
shown. Curves are fits to the power law function. Hash marks 
at the top indicate bin boundaries for > 1.5 GeV/c. 



III. RESULTS 



ofNbin = B-Np^^t, were obtained by fit to be 1.41±0.08, 
1.34±0.08, and 1.38±0.08 for HIJING, VENUS, and MC 
Glauber. The scaling exponents a for these models are 
approximately 4/3 due to the fact that Npart oc and 

It is worthwhile to note here that distribution differ- 
ences among HIJING and VENUS are mainly due to dif- 
ferent overlap functions. The Npart and Nhm distribu- 
tions from them are nearly identical if the same overlap 
functions are used in these two model calculations. Fig.O 
shows that over a broad range of centrality the model de- 
pendent uncertainties of Npart and Ntin are within 10% 
and 20%, respectively. 



Inclusive pt distributions of (/i+ + h~)/2 within |?7| < 
0.5 have been published previously Q. The independent 
analysis reported in this article shows that the differences 
to the published results for all measured px points are 
within 10%, which is comparable to the systematic uncer- 
tainties for Pt < 2 GeV/c and is less than the systematic 
uncertainties for the high pT region. Fig. 0] shows inclu- 
sive Pt distributions of {h~^ + h~)/2 within 0.5 < j??! < 1 
for various centrality bins. The error bars are the quadra- 
ture sum of statistical error and systematic uncertainty, 
and are dominated by the latter except for the highest 
Pt point in the peripheral bins. The curves in Fig.0]are 
power law function (Eq. ^ fits to the spectra. 
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FIG. 5: Ratios of pr distributions within 0.5 < |r;| < 1 to 
those within |?7| < 0.5 in various centrality bins. Error bars 
show statistical errors while caps are the quadrature sum of 
statistical errors and remaining systematic uncertainties. 



Fig. 13 shows ratios of px distributions within 0.5 < 
\r]\ < 1 to those within < 0.5 in various centrality 
bins. Note that Fig. Inland the succeeding figures utihze 
the Pt distributions within \r]\ < 0.5 obtained here. Us- 
ing identical cuts and correction procedures across the 
full pseudorapidity region minimizes the systematic un- 
certainties in the relative comparisons. The error bars 
in Fig. |31 show statistical errors only while the caps are 
the quadrature sum of statistical errors and systematic 
uncertainties which cannot be canceled out. Remaining 
systematic uncertainty includes the variation due to track 
quality cuts, the uncertainties of the px smearing correc- 
tions for the two 77 regions, and the partial uncertainty of 
background subtraction related to the ?7-dependcnt part 
discussed in Sec. II. 

Fig. El shows the same ratio of 0.5 < |ry| < 1 to 
\ri\ < 0.5 in the 0-5% most central bin. The points are 
our measurements and the error bars include statistical 
and remaining systematic uncertainties. The solid curve 
is the same ratio from PYTHIA calculations for 130 
GeV p + p collisions. Other curves are ratios from HI- 
JING predictions of 130 GeV Au -I- Au collisions without 
shadowing and without quenching (dotted curve), with 
shadowing and without quenching (dashed curve), and 
with shadowing and with partonic energy loss being 2.0 
GeV/fm (dot-dashed curve). The results show that the 
effects on the pseudorapidity dependence of both nuclear 
shadowing and partonic energy loss as implemented in 
HIJING are too small to be tested in the measured kine- 
matic region under current experimental uncertainties. 

No significant differences are observed in the compar- 
isons of the inclusive charged hadron yields between the 
two T] regions in Fig. [SI and Fig. over a broad range of 
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FIG. 6: Ratio of pr distribution within 0.5 < < 1 to that 
within |r;| < 0.5 in the 0-5% most central bin. Points are 
measurements and error bars include statistical and remaining 
systematic uncertainties. Curves are described in the text. 
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FIG. 7; dN/dr] distributions for pr > 2 GeV/c and -1 < 
77 < 1 scaled by Nbin and divided by the A'^iV reference. 



centrality for all measured px points. It suggests that an 
approximate boost invariant condition might be estab- 
lished in the early stage of collisions. The suppression 
pattern of the particle yield has little 77 dependence in 
the measured region though the particle yield itself is 
sensitive to partonic energy loss. A measurement of this 
ratio between 77 = 2.2 and = from the BRAHMS Col- 
laboration shows that the ratio is below unity at ^ 4 
GeV/c 113. 

Fig. [7| shows dN/df] distributions for px > 2 GeV/c 
and — 1 < 77 < 1 in various centrality bins. The er- 
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ror bars show statistical uncertainties while the caps are 
the quadrature sum of statistical and systematic uncer- 
tainties. The systematic uncertainties are dominant and 
highly correlated. The dN/drj distributions are scaled by 
Nbin and divided by the NN reference. Due to nearly 
complete 77-independence of the NN reference data for 
Pt > 0.2 GeV/c within — 1 < ry < 1 as shown in Fig. [21 
and Fig. El constant dN/ di] of the NN reference are used 
in Fig. Therefore, the shapes of the dN/dr] distri- 
butions for the Au -I- Au collisions are preserved. The 
uncertainties on both N^m and the NN reference data 
are shown in the shaded regions around the lines at unity 
which represent the binary collision scaling. Ratios below 
unity in the figure show that the high px hadrons over 
2 GeV/c are suppressed with respect to those in p + p 
collisions. The shape of dN/drj for the high px hadrons 
is nearly flat. No significant centrality dependence of the 
dN/drj shapes within — 1 < 77 < 1 is observed. Similar 
behaviors are observed for px > ^ GeV/c except larger 
suppressions in the central bins. For example, in the 0- 
5% most central bin, the average ratio is 0.41 ± 0.10 for 
PT > 4 GeV/c while it is 0.64 ± 0.10 for pr > 2 GeV/c. 



IV. DISCUSSION 
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FIG. 8; Upper panel: ratio of charged hadron yields within 
I77I < 1 for Au + Au relative to the A'^A'' reference, scaled by 
Npart/'i as a function of centrality for a pr bin at pr = 2.05 
GeV/c. The curve is fit to -B- A^a^j. Lower panel: participant 
scaling exponent v of charged hadron yields as a function of 
Pt within \rj\ < 1. 



The charged hadron yield per participant pair at 
y/Sj,, „ = 130 GeV shows a slow increase as a function 



of iV, 



part 



35l l36l . Such slow increase of hadron multi- 



plicity as a function of centrality at RHIC has been con- 
sidered by Kharzeev et al. [3^ in the framework of parton 
saturation. They argued that the hadron multiplicity as 
a function of centrality would increase faster if produced 
jets lose energy radiating soft gluons that in turn frag- 
ment into hadrons at mid-rapidity. As a result of the 
parton saturation, it is predicted that hadron multiplic- 
ity should scale with Npart at a moderately high px (up 
to 6-8 GeV/c at RHIC energies). An explanation of the 
slower than expected increase in fragmentation models is 
that the effective energy loss is significantly reduced in a 
thermal environment due to detailed balance |^. Recent 
experimental results in c? -I- Au collisions at y/s^ = 200 
GeV support the idea that the suppression of high px 
hadron production in Au -I- Au collisions at mid-rapidity 
is due to final state interactions rather than parton sat- 
uration in the initial state Hj, jUl . 



In the upper panel of Fig. |S1 we plot charged hadron 
yield per participant pair within \rj\ < 1 normalized to 
that of nucleon-nucleon collisions as a function of Npart 
for Px = 2.05 GeV/c. The error bars are the uncertain- 
ties of data while the caps are the quadrature sum of the 
uncertainties of both data and Npart- The shaded regions 
around unity show systematic uncertainties of the NN 
reference data. The result shows that the ratio is above 
unity and increases with Npart- 

Dependence of the charged hadron yield on Npart can 



be studied by fitting the yield by the following function 
d^N 



dpxdij 



B ■ Npart 



(4) 



in different px bins. Such an example is shown as a curve 
in the upper panel of Fig. [Slfor px = 2.05 GeV/c. The 
fit parameter, v{px), is given in the lower panel of Fig. (HI 
as a function of px- The error bars are the uncertainties 
of the fit parameters associated with the uncertainties 
of data. The lines and shaded regions are binary colli- 
sion (Nbin) and participant {Npart) scaling exponents and 
uncertainties to Npart- No clear evidence of participant 
scaling over the whole measured px region is observed. 
The approximate participant scaling of the hadron yield 
at high Px observed by PHOBOS 39] appears to be ac- 
cidental. 

In a scenario with continuous energy loss of particles 
through a medium, the energy loss would lead to a shift 
in the mean px of these particles. If the energy loss 
contributes to additional particle production in the low 
Px region, the mean px of low px particles would also be 
modified. The truncated mean px, defined as 



iPx^niPT') 



/pTut Pt ■ dN/ dpx 



dpi 



dN/dpx 



dpi 



cut 

Px 1 



(5) 



is used to study the variation of mean px as a function 
of Px scale with respect to NN reference data. Fig. (21 
shows the truncated mean px ratios between Au -I- Au 
and p + p collisions as a function of p™* for central (0- 
5%) and peripheral (60-80%) collisions. The errors are 
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binary scalings. In two-component particle production 
models, the binary scaling has been associated with hard 
parton scatterings and the participant scaling with the 
soft processes. In our study we empirically decompose 
the particle yield into N^art and Nun scaling compo- 
nents, i.e. 



N, 



part 



(6) 



where npp{p^*) and are the hadron multiplicity 

and the fraction of particle yield attributable to hard 
processes in p + p collisions, respectively. Fig. 1101 shows 
the binary scaling fraction in Au -I- Au collisions, defined 
as 



Fip: 



cut -J 



dN/drjipf^) 



(7) 



Note that F{p^*) does not depend on npp(p™*) since 
both numerator and denominator of Eq. [7| contain 
npp{p^*). There is a distinguishable trend as a func- 
tion of Npart from p^* = 3.8 GeV/c to lower pf^. This 
trend is consistent with the i'{pt) dependence in Fig.|Sl 
It is worth noting that F m 70% in central Au -I- Au 
collisions at p™* = 3.8 GeV/c. However, one should ex- 
ercise caution when relating this fraction to hard parton 
scattering processes, particularly at lower px where high 
Pt particles may suffer large energy losses in the medium 
and become soft. 



0' ' ' 

100 200 300 100 200 300 

part 

FIG. 10: Binary scaling fraction in pt > p'r^ within j?;} < 1 
as a function of centrality for selected p^"*. For p™* > 2 
GeV/c, the fraction F decreases with centrality. 

combined statistical and systematic uncertainties while 
the caps are the quadrature sum of the uncertainties of 
both the Au -f Au data and the reference data. 

In peripheral collisions at high pT (p^"* > 3 GeV/c) 
the truncated mean pT of particles is approximately the 
same as for p -|- p collisions above the same p™* (Fig. EJl . 
The ratio in the low pT region is above unity indicating 
the effects of the Cronin effect and/or radial flow in pe- 
ripheral collisions. For central collisions, the truncated 
mean pT for p™* > 3 GeV/c is approximately 15% lower 
than the truncated mean pT from p -V p collisions at the 
same p™*, consistent with the scenario for partonic en- 
ergy loss in this pt region. The significantly larger ratio 
in the low pt region probably reflects the combined ef- 
fects of larger radial flow, the Cronin effect, and pt shift 
of particles due to energy loss, which cannot be decoupled 
with the present data. 

Fig. |S1 (bottom panel) indicates that over a broad pt 
region particle production falls between participant and 



V. CONCLUSION 

We have presented inclusive distributions of (Ji^ -t- 
/i^)/2 from STAR at RHIC in the region 0.5 <\t]\<\ 
and compared them to distributions for |7/| < 0.5, finding 
no significant differences in the region of 0.2 < < 6.0 
GeV/c. We find that the dN/drj distributions for —1 < 
77 < 1 are nearly flat for all centralities. The charged 
hadron yield as a function of px shows no clear par- 
ticipant scaling in the measured pt region. The binary 
scaling fraction in the two-component model shows a de- 
crease with centrality for p™* > 2 GeV/c and is about 
70% at p"^* = 3.8 GeV/c for central colhsions. 
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